PARTIAL REGULARITY OF SOLUTIONS TO THE 
FOUR-DIMENSIONAL NAVIER-STOKES EQUATIONS 



HONGJIE DONG AND XUMIN GU 

Abstract. In this paper, we consider suitable weak solutions of incom- 
pressible Navier-Stokes equations in four spatial dimensions. We prove 
that the two-dimensional time-space Hausdorff measure of the set of 
singular points is equal to zero. 



1. Introduction 

In this paper we consider the incompressible Navier-Stokes equations in 
four spatial dimensions with unit viscosity and an external force: 

ut + u-Vu- Au + Vp = f (1.1) 
V-n = (1.2) 

in a bounded cylindrical domain Qx = ^2 x (0,T), where Cl C M^. We 
are interested in the partial regularity of suitable weak solutions {u,p) to 

(HID-dEl]). 

We say that a pair of functions {u,p) is a suitable weak solution to 
([riI-([L2]) in Qr if n e LooiO,T; L2{n;R'^)) n L2{0,T-W^{n;R^)) and p G 
-^3/2 (Qt) satisfy (|l.ip - (|1.2p in the weak sense and additionally the gen- 
eralized local energy inequality holds for all non-negative functions ip S 



ess sup / |n(x,s)| ^{x,s)dx + 2 / |Vu| ipdxds 
o<s<t Jn JQt 

< / \u\'^{ipt + Aip) + {\u\'^ + 2p)u-V7p + f-mpdxds. (1.3) 

We will prove that for any suitable weak solution {u,p), the two dimen- 
sional space-time Hausdorff measure of the set of singular points is equal to 
zero. 

The problem of the global regularity of solutions to the Navier-Stokes 
equations in three and higher space dimensions is a fundamental question 
in fluid dynamics and is still widely open. Meanwhile, many authors have 
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studied the partial regularity of solutions. In three dimensional case, Schef- 
fer established various results for weak solutions in [12^ [T3| . In a celebrated 
paper [Tj, Caffarelli, Kohn, and Nirenberg firstly introduced the notion of 
suitable weak solutions, which satisfy a local energy inequality. They proved 
that for any suitable weak solution, there is an open subset where the ve- 
locity field u is regular and the ID Hausdorff measure of the complement of 
this subset is equal to zero. In [TT], Lin gave a more direct and simplified 
proof of Caffarelli, Kohn and Nirenberg's result with zero external force. La- 
dyzhenskaya and Seregin gave a detailed account in [9] later. We also refer 
the reader to Tian and Xin |17) . He [5], Seregin [15], Gustafson, Kang, and 
Tsai [1], Vasseur [18J, Kukavica [8], and the references therein for extended 
results. 

For the four or higher dimensional Navier-Stokes equations, the problem 
is more super-critical. In [141, Scheffer showed that there exists a weak 
solution u in X M+, which may not satisfy the local energy inequality, 
such that u is continuous except for a set whose 3D Hausdorff measure 
is finite. In [2], Dong and Du showed that for any local-in-time smooth 
solution to 4D Navier-Stokes equations, the 2D Hausdorff measure of the 
set of singular points at the first potential blow-up time is equal to zero. 
Moreover, for stationary high dimensional Navier-Stokes equations, Struwe 
[16j proved that suitable weak solutions are regular outside a singular set 
of zero ID Hausdorff measure in R^, and Kang [6j improved Struwe's result 
up to the boundary for a smooth domain C M^. Recently, Dong and 
Strain [3] studied the partial regularity for suitable weak solutions of 6D 
stationary Navier-Stokes equations, and proved that solutions are regular 
outside a singular set of zero 2D Hausdorff measure. Based on Campanato's 
approach, the main idea in [3] is to first establish a weak decay estimate 
of certain scaling invariant quantities, and then successively improve this 
decay estimate by a bootstrap argument and the elliptic regularity theory. 

Because time corresponds two space dimensions, in some sense the 4D 
non-stationary Navier-Stokes equations is similar to 6D stationary prob- 
lem. Given the result in [3], it is natural to ask whether Caffarelli-Kohn- 
Nirenberg's theorem can be extended the 4D non-stationary case. Here the 
main difficulty stems from the fact that certain compactness arguments ap- 
peared, for instance, in the original paper [l] as well as [HI [9] break down 
in the 4D case. We note that the results obtained in the [2] cannot be con- 
sidered as a genuine extension of the theorem, as the set of singular points 
is only estimated at the first blow-up time for local smooth solutions. The 
objective of this paper is to give a complete answer to this question. 

We state our main results, where we use some notation introduced at the 
beginning of the next section. 

Theorem 1.1. Let Q be a open set in W^, f G Lq ioc{Qt)- Let {u,p) be 
a suitable weak solution of (|l.ip - (jl.2p in Qt- There is a positive number 
Eq satisfying the following property. Assume that for point a zq G Qt, the 
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inequality 

limsup E(r) < Eq 
holds. Then zq is a regular point. 

Theorem 1.2. Let Q be a open set m / G Lq^i^(.{Qt)- Let {u,p) be a 
suitable weak solution of (|1.1|) - (|1.2|) in Qj-. There is a positive number Eq 
satisfying the following property. Assume that for a point zq £ Qt and for 
some po > such that Q{zq,pq) C Qt and 

C{po) + D{po) + F{po) + G{po) < so. 

Then zq is a regular point. 

Theorem 1.3. Let Q be a open set in M^, / S LqJ(,c{Qt)- Let {u,p) be a 
suitable weak solution of (jl.ip - (jl.2p in Qt. Then the 2D Hausdorff measure 
of the set of singular points in Qt is equal to zero. 

Compared to ^ , as mentioned above Dong and Du concerned with local- 
in-time smooth solution to 4D Navier-Stokes equations with zero external 
force. Due to the lack of compactness, they used Schoen's trick in their proof. 
In our paper, we shall consider suitable weak solutions, and thus Schoen's 
trick is no longer applicable. Our proofs exploit the aforementioned idea in 
[3] and use Campanato's approach. There are two main differences between 
our problem with the one in [3]. The first one is that we do not have the 
same end-point Sobolev embedding inequality which was used in [3j. To this 
end, we introduce an additional scale-invariant quantity F, which is a mixed 
space-time norm of the pressure p, and use an interpolation inequality. As 
a consequence, we cannot archive the same optimal decay rate as in [3]. 
Nevertheless, it turns out that the decay rate, although not optimal, still 
suffices for our purpose in the subsequent step. The other difference is 
that, as our problem is time-dependent, we cannot use the elliptic regularity 
theory to improve the decay rate in the final step as in [3]. Naturally, we 
appeal to the parabolic regularity theory instead as well as a Poincare type 
inequality for solutions to divergence form parabolic equations. 

We remark that by using the same method we can get an alternative proof 
of Caffarelli-Kohn-Nirenberg's theorem for the 3D Navier-Stokes equations 
without using any compactness argument. It remains an interesting open 
problem whether a similar result can be obtained for five or higher dimen- 
sional non-stationary Navier-Stokes equations. It seems to us that four is 
the highest dimension to which our approach (or any existing approach) 
applies. In fact, by the imbedding theorem, we have 

Loo((0, T); n L2((0, T); ^ L2(rf+2)/d((0, T) x J]), 

which implies nonlinear term in the energy inequality cannot be controlled 
by the energy norm alone when d > 5. 

We organize this paper as follows: In Section [21 we introduce the notions 
of some scaling invariant quantities and essential settings which would be 
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used throughout the paper. In Section [3l we prove our results in three steps. 
In the first step, we give some estimates of the scaling invariant quantities, 
which essentially follows the argument in [2]. In the second step, we establish 
a weak decay estimate of certain scaling invariant quantities based on the 
estimate we proved in the first step by using an iteration metho. In the last 
step, we improve the decay estimate by a bootstrap argument, and apply 
parabolic regularity to get a good estimate of the L3/2-mean oscillations 
of u, which yields the Holder continuity of u according to Campanato's 
characterization of Holder continuous functions. 



In this section, we will introduce the notation which would be used 
throughout the article. Let be a domain in some finite-dimensional space. 
Denote Lp(J^;R") and Wp^(J^;M") to be the usual Lebesgue and Sobolev 
spaces of functions from into M". Let p £ (1, oo) and — oo<S'<r<oo. 
We denote Tip to be the solution spaces for divergence form parabolic equa- 
tions. Precisely, 



nlin X {S,T)) ={u:u,Du£ Lp{n X {S,T)), ut G W^^i^ x (5,r)}, 



inf {11/11 Lj,(nx(S,T)) + I|5||Lp{Qx(5,T)) I V = V • 5 + /} < oo. 

We shall use the following notation of spheres, balls, parabolic cylinders, 
and parabolic boundary 

5(xo,r) = {xGM^||x-xo| =r}, 5(r) = 5(0,r), 5 = 5(1); 
^(xcr) = {xeM^||x-Xo| <r}, B(r)=5(0,r), B = B{1); 
Q{zo, r) = B{xo, r) x {to - r^, to), Q{r) = Q(0, r), Q = Q(l); 
dQizo,r) = {S{xo,r) x [to - ,to)) U {{to - ,y) \y G B{xo,r)}, 

where zo = (xq, to)- 

We also denote mean values of summable functions as follows: 



2. Notation and Settings 



where EIp^(ri x (5, T)) is the space consisting of all generalized functions v 
satisfying 




Here 1^41 as usual denotes the Lebesgue measure of the set A. 
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Now, we introduce the following quantities: 

|2 



A{r 
E{r 
C{r 
D{r 
F{r 
G{r 



1 



1 



A{r, zq 
E{r, Zq 
C{r, Zq 

D{r,zo 
F{r,zo 
G{r,zo) = 



1 

ess sup — 

to-r^<t<to JB{xo,r) 



\u{x,t)\ dx, 



\Vu\^dz, 



\u\^ dz, 



\p ■ 



1 



f 

Q{zo,t) 



\p\xo,r\ ^ dz, 



B{xo,r) 



\p-\pUA^^'"dx) 2" dt 



2a 
1 + a 



\frdz 



Q{zo,r) 



where a G (0, 1) is a number to be specified later. Notice that all these 
quantities are invariant under the natural scaling: 

ux{x,t) = Xu{Xx,X^t), px{x,t) = \^p{\x,\^t), fxix,t) = A^/(Ax,A^t). 

We are going to estimate them in Section [3l Note that the quantity F is 
auxiliary and will only be used in the first two steps of the proof in order to 
give a weak decay estimate of other quantities. 

We finish this short section by introducing a pressure decomposition which 
would play a important role in our proof. Let r]{x) be a smooth function on 
supported in the unit ball B{1), < rj < 1 and rj = 1 on B{2/3). Let zq 
be a given point in Qt and r > a real number such that Q(zQ,r) C Qt- 
It's known that for a.e. t G (to — r'^,'to), in the sense of distribution, one has 



Ap 



dxidxj 



{uiUj) + V • / 



52 



This will hold for a weak solution to (|l.ip - (jl.2p . For these t, we consider the 
decomposition 

P = Pxo,r + hxo,r in B{xo,r), (2.1) 
where Pxo,r is the Newtonian potential of 

g^.g^ , ((^i-K]xo,r/2)(^i-K-]xo,r/2)??(2(x-a:o)/r))+V-(/?7(2(x-xo)/r)). 



Then /ixo,r is harmonic in B(xo,r /3). 
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3. The proof 

In our proof of the results, we will make use of the following well-known 
interpolation inequality. 

Lemma 3.1. For any function u E VF2^(M^) and real numbers q £ [2,4] and 
r > 0, 

\u\'' dx < N{q) [ I \Vu\'^dxy~^{ f {ul"^ dx)' 

J Bj- J Bf 

+ r-2('?~2)( I \u\^dxy'^ 

J Br 

Let C := dt — dx^iaijdx^) be a (possibly degenerate) divergence form 
parabolic operator with measurable coefficients which are bounded by a 
constant K > Q. We will use the following Poincare type inequality for 
solutions to parabolic equations. See, for instance, [2, Lemma 3.1]. 

Lemma 3.2. Let zq G M'^+^ p G (l,oo), r G (0,oo), u G V-l^iod^'^'"^) , 
g = {gi, . . . , Qd), f G Lpjoci^'^^^)- Suppose that u is a weak solution to 
Cu = V ■ g + f in Q{zo,r). Then we have 



\u{t, x) - {u),,^Ry dz < NrP / (I Vn|P + + rP|/|P) dz, 

where N = N{d,K,p). 

Now we prove the main theorems in three steps. 

3.1. Step 1. First, we control the quantities A,C,D,F in a smaller ball 
by their values in a larger ball under the assumption that E is sufficiently 
small. Here we follow the arguments in [2], which in turn used some ideas 
in pin]. 

Lemma 3.3. Suppose 7 G (0,1), p > are constants and Q{zq, p) C Qt- 
Then we have 

Chp) < N[^-'A'/\p)E{p) + ^-'/'A'/Hp)E'/\p) + 7C(p)], (3.1) 

where N is a constant independent 0/7,^ and zq. 

The proof can be found in [2]. 

Lemma 3.4. Suppose a G [1/11, 1/2], 7 G (0, l/8],/9 > are constants and 
Q{zq,p) C Qt- Then for any z\ G Q{zQ^p/S) we have 

F{^p,z^)<N{a)[r^ATTt{^^)E^{^)+^'^F{p) + r^G^/\p)l (3.2) 

where N{a) is a constant independent of 7,p and zq. In particular, for 
a = 1/2 we have 

Di^p, zi) < iV[7-3AV2(^)i^(|) + ^5/2 D{p) + ^-^G^/\p)]. (3.3) 
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Moreover, it holds that 



(3.4) 

Proof. First, we assume 1/3 < a < 1/2. Denote r = jp. Recall the de- 
composition of p introduced in (12. ip and the definition of rj. By using the 
Calderon-Zygmund estimate, Lemma l3.ll with q = 2(1 + a), the Poincare 
inequality, and the Sobolev embedding inequality, one has 



B{xo,r) 

<N [ \u- M,„,,/2p(i+") dx + N [ |A-iV • (/r/(2(x - xo)/r))| 

JB{xo,r/2) JB{xo,r) 



1-a 



<N{[ iVul^dxf'^il |u-H,„,,/2pdx) 

JB{xo,r/2) JB{xo,r/2) 

+ iVr-4"( / |^_[^^^^/2l'dx)^+" + iV( / Ifl^dx)"^ 

Jb(xo,t/2) JB{xo,t/2) 



< N( 



B{xo,r/2) 



|Vn|2dx)'"( / \u\^dxy-'' + N{ [ \f\^dx) 

JB{xo,r/2) JB{xo,r/2) 

(3.5) 



where > 1- Here we also used the obvious inequality 

/ \U- W\xo,r/2?dx < I 

JB(xn.r/2) JE 

Similarly, 



|up dx. 

B{xo,r/2) ' ' JB{xo,r/2) 



I \Pxo,p\^-^''dx<N{ [ \Vu\^dx)^''{ [ |n|2dx)'"" 

Jb{xo,p) Jb{xo,p/2) Jb{xo,p/2) 



+ iV( / |/|Wdx) 4 . (3.6) 

Jb{xo,p/2) 



Since h^^^p is harmonic in B(xo, p/3), any Sobolev norm of /ixo,p iii a 
smaller ball can be estimated by its Lp norm in B(xo, p/3) for any p £ [1, oo]. 
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Thus, by using the Poincare inequahty one can obtain, for a.e. t, 
/ \h -\h ] dr 

I \'''Xo,p ["'Xo,p\xi,r\ ""^ 

JB(xi,r) 

J B{xi,t) 

Ih 

'XO,p\ 



<^^5+a |y, |l+a 



B(xi,r) 
P Jb{xo,p/3) 

< iV(-)'+"[ / \pix,t) - [p].o,pl'+° + |Pxo,p(x,t)|^+"dx]. (3.7) 

P Jb{xo,p) 

Since Pxo,r + Kq^t = p = Pxo,p + Ko,p in B{xi,r), from ([33]), ([M]) and 
p.7p we get, for a.e. t. 



J_B(a;i,r) 



' dx 



— / |P2;o,p [Pxo,p]a;i,r I dx -\- I \hxo,p [^a::o,p]a::i,r| dx 
JB{xi,r) JB{xi,r) 

< I \f) \^^"- dr 4- / \h —\h\ \^^'^ dr 

— I \l'xo,p\ "'■^ / \"'Xo,p ["'Xo,p\xi,r\ "'■^ 

JB{xi,r) Jb(xi,t) 

<N{[ |Vt/pdx)'"( / lul'dxf-" 

Jb(xo,p/2) Jb{xq,p/2) 

+ iV(-)'+" / \p{x,t)-\pUJ^+^dx + N{ [ l/lWdx)"^. 

P JB(xQ,p) JB{xo,p) 

(3.8) 

Raising to power — , integrating with respect to t in (ti — r^, ti), and using 
2a 



Holder's inequahty complete the proof of (|3.2p and (|3.3p . 

Now in the case 1/11 < a < 1/3, we cannot use the Sobolev embedding 
inequality directly in (j3.5p because < 1. However, since function rj has 
compact support, by using Holder inequality, we can get 

nI IA-^V • (/r/(2(x -xo)/r))|^+"dx 

JB{xQ,r) 

<Ar/(i-w)( I |A-iV-(/r?(2(x-xo)/r))|^+'^dx)^ 



B(xo,r) 

A(^ i+^V f . ,±tM 

<Nr^^~~'{ \f\~ dx) 4 

JB{xQ,r/2) 



where 1/3 < /3 < 1/2. Noting that < g and ^ < 6, we then prove 

(j3.2p in the same way as the case 1/3 < a < 1/2. 
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To prove ()3.4p . we use a slightly different estimate from ()3.7p . Again, 
since h is harmonic in B{xq, p/3), we have 

h _ f/l 1 |3/2 fjrf. 

"'Xo,p l"'Xo,plxi,r\ "'■^ 

B{xi,r) 

JB{x-i,r) 
B{xur) 



^ ^ 3/2+6/(l+») [ / l^-0,p(^,0 - b]xo,pl'+"]^ 



^11/2 



JB{a;o,p) 

Similar to p.Sh . we obtain, for a.e. t, 

B{xi,r) 



<N{ [ \Vu\'^dx){ [ \u\'^dx)^^^ 

Jb{xo,p/2) Jb{xo,p/2) 
K L J B(xo,p) 



2(l-t-a) 



/■ 3a /■ 3(l-a) 

+ (/ |Vu(x,t)|2(i3;)i+"( / |M(x,t)p(ix)2(i+«) 

J B{xo,p) J B{xQ,p) 



yB(xo,p) 

Integrating with respect to t in {ti — r^,ti) and applying Holder inequality 
completes the proof of p.4j) . □ 

Lemma 3.5. Suppose 6 G (0, 1/2], p > are constants and Q{zo,p) C Qt- 
T/ien we have 

A{ep) + < NO-^[C^/\p) + + C^/^ip)D^/\p) + Gip)] . 

In particular, when 6 = 1/2 we have 

A{p/2) + E{p/2) < N[C'/'{p) + C{p) + C^''{p)D^/\p) + G{p)\ . (3.9) 

Proof. Let r = Op. In the energy inequality p.3p . we set i = to and choose 
a suitable smooth cut-off function -0 such that 

= in QtQ\Q{zo,p), < -0 < 1 in Qt, 

xP = l in Q{zo,r), \V^P\ < Np-\ \dtiP\ + \V^iP\ < Np-^ in Qt,. 
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By using (ll.Sp and because u is divergence free, we get 



N [ 1 f If 
A{r) + 2E{r)<^ ^ \u\Uz + - i\u\^ + 2\p-[pl,J)\u\dz 

Lp^ jQ(^„,p) P Jq(zo,p) 



+ / \fMdz 

IQ(zo,p) 



Using Holder's inequality and Young's inequality, one can obtain 

/ \u\Uz < ( / \ufdzf'{ f dzY^' < p^C^/\p), 

jQ(zo,p) JQ{zo,p) Jq{z(),p) 

[ \p-[pUJ\u\dz < ( [ \p-[pUjy'dzf\ I \u\^dzf'^ 
Jq{zo,p) Jq(zo,p) Jq(zo,p) 



and 



/ \f\\u\dz<p^ [ \f\^dz + ^ [ 

JQ(zn,p) JQ(zn,p) P JQ\ 



luP dz 



'Q{zo,p) JQ{z„,p) P jQ(zo,p) 

<p'\[ \ffdz 



+ ^ \u\ dz. 



P'^ Jq{zo,p) 



Then the conclusion follows immediately. 



As a conclusion, we can obtain 



□ 



Proposition 3.6. For any eq > 0, there exists ei > small such that for 
any zq £ Qt satisfying 

limsup£;(r) < £1, (3.10) 

we can find po sufficiently small such that 

A{po) + Eipo) + C(po) + D{po) + F{po) < sq. (3.11) 

Proof. First, we prove (j3.1ip without the presence of F on the left-hand 
side. For a given point 

ZQ = {xo,to) £ Qt 

satisfying (|3.10p . choose pi > such that Q{zo,pi) C Qt- Then for any 
p £ {0,pi] and 7 G (0, 1/8), by using ()3.9p and Young's inequality, 

Ai^p) + E{^p) < N[C^/\2^p) + C(27p) + D{2^p) + G{2-fp)\ . 
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Then, combining with (I3.ip and (I3.3p and using Young's inequahty again, 
we have 

A{-fp) + E{-ip) + C(7p) + D{-ip) 

< iV[72/3c2/3(p) + f + ^C{p) + -fA{p)] 

+ iV7-ioo(S(p) + E^ip) + G{p)) + N^^/^ 

< N^^/^ [Aip] + Eip) + C{p) + D{p)] + N^^/^ 

+ Nj-'^'{E{p) + eHp) + G{p)). (3.12) 
Since / G -^^6,1oc(Qt); we have 

G{P) < ll/llL(g(.o,p,))P^ (3.13) 

It is easy to see that for any eq > 0, there are sufficiently smah real numbers 
7 < l/(2A^)3/2 and Ei such that if ()3.10p holds then for all small p we have 

Ar^2/3 + Nj-''^{E{p) + E-'ip) + G{p)) < eo/2. 

By using (j3.12p . we can obtain 

A{po) + C{po) + D{po) < eo 

for some po > small enough. To include F in the estimate, it suffices to 
use (1321). □ 



3.2. Step 2. In the second step, first we will estimate the values of A, E, 
C, and F in a smaller ball by the values of themselves in a larger ball. 

Lemma 3.7. Suppose p > 0, 9 £ {0, 1/16] are constants and Q{zi,p) C Qt- 
Then we have 

A{9p) + E{9p) < N9^A{p) + NQ-"^ [A{p) + E{p) + F{p)] ^'^ 

+ N9-^G^/\p) + Ne-^G{p), (3.14) 
where N is a constant independent of p, 9, and zi. 
Proof. Let r = 9p. Define the backward heat kernel as 

1 



^ ' 47r2(r2 + ti -t)2 



In the energy inequality (jl.3p we put t = ti and choose = r^/" '■= 
T(j)i[x)(t)2{t), where 0i, 4>2 are suitable smooth cut-off functions satisfying 

(/)i = in M.\B{xi,p), 0<(/)i<l in M^, (l)i = l in B{xi,p/2), 

(t>2 = in (-oo,ti - p^) U (ti +/)^,+oo), < (/>2 < 1 in M, 

^2 = 1 in (ti - /)V4,ti +pV4), \dtH < N p-^ in M, 

|V0i| < 7V/3-Mv2,/.i| < in M^. 

(3.15) 
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By using the equality 

Ar + Ft = 0, 

we have 

\u{x,t)\'^T{t,x)(t>{x,t)dx + 2 / \\7u\^T4>dz 
B{xup) Jq{zi,p) 

< [ {\u\'^iT(j)t + TA(j) + 2V(l)Vr) + {\u\'^ + 2p)u-iTV<J) + (j)VT)}dz 



+ / \f\\u\\r<p\dz. (3.16) 

With straightforward computations, it is easy to see the following three 
properties: (i) For some constant c > 0, on Q{zi,r) it holds that 

T(j) = T> cr""^. 

(ii) For any z £ Q{zi,p), we have 

\r{z)(l){z)\ < Nr~^, |0(z)VF(z)| + |V0(z)F(z)| < Nr'^. 

(iii) For any z £ Q{zi, p)\Q{zi,r), we have 

|F(z)</)t(z)| + |F(z)A</)(z)| + |V0VF| < TVp-^. 
These properties together with (j3.15p and (|3.16p yield 

A{r) + E{r) < N[e^A{p) + e-^{C{p) + D{p)) + e~^G{p)\ . (3.17) 
Owing to (|3.1|) with g = 3, we can get 

C{p/S) < NC{p) < N[A{p) + E{p)f\ (3.18) 
By using (|3.4p with 7 = 1/8, we have 

Dip/8) < N[A{p) + Eip) + F{p)f' + NG'/\p). (3.19) 

Upon combining (|3.17p (with p/8 in place of p) to (|3.19p together, the lemma 
is proved. □ 

Lemma 3.8. Suppose p > is constant and Q{zi,p) C Qt- Then we can 
find 9i G (0, 1/256] small such that 

Aie.p) + E{e,p) + F{9,p) < i [Aip) + Eip) + F{p)] 

+ N{0^) [A{p) + E{p) + F{p)] ^'^ + N{e,) [G{p) + G^l\p)\ , (3.20) 
where N is a constant independent of p and zi . 
Proof. Due to ([32]) and (l3Jl]l . for any 7,6* G (0, 1/16], we have 

Fi-fOp) <N[j-\A{ep) + E{ep)) + 7(3-")/(i+")F(0p) + 7-2^1/2(0^)] 
<iV7'2^2^(/>) +7(3-")/{i+")0-2f(p) +7-2g1/2(^) 

+ Nj-'e-'[Aip) + E{p) + F{p)f\ (3.21) 
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and 

A{^ep) + E{-fdp) < N{jefA{p) + Ni-fO)-^ [A{p) + E{p) + F{p)] 

+ Ni-fOy^Gip) + iV(7e)-3G3/4(p). (3.22) 

Now we set a = 1/5 such that (3 — a)/(l + a) = 7/3 > 2, and choose and 
fix 7 and sufficiently smah such that 

Upon adding (|3.21|) and (|3.22|) . we can obtain 

Aijep) + Ei^ep) + Fi^ep) < ^A^p) + n [a{p) + e{p) + f{p)] 

+ N[G{p) + G'l\p) + G'/\p)l 

where only depends on 6 and 7. After putting 61 = jO, the lemma is 
proved. □ 

In the next proposition we will study the decay property of A, C,E and 
F as the radius p goes to zero. 

Proposition 3.9. There exists Eq > satisfying the following property. 
Suppose that for some zq £ Qt and po > satisfying Q{zo,po) C Qt we 
have 

C(po) + Dipo) + F{po) + G(po) < £0- (3.23) 

Then we can find N > and oq £ (0, 1) such that for any p £ (0, po/8) and 
z\ G Q{zq, p/S), the following inequality will hold uniformly 

A{p, zi) + G{p, zi) + Eip, zi) + F{p, zi) + D{p, zi) < N p^\ (3.24) 

where N is a positive constant independent of p and z\ . 

Proof. Fix the constant 9i £ (0,1/256] from Lemma 13.81 and let N{6i) > 
is the same constant from (j3.20p . Due to (j3.9p . ()3.23p . and (|3.2p . we first 
choose e' > and then eo = ^oi^') > sufficiently small such that, 

N{9i)^ < 1/4, iV(0i) (eo + el^^) < e' /2, (3.25) 

A{po/2) + E{po/2) < ^, 
and for any zi £ Q{zq, pq/8), 

Fipo/8,z,)<N[Ay^P^)E'/H'-^) + F{po) + G'^\po)] < 
By using 

Q{zi,po/8) C Q{zo,po/2) c Qt, 

we then have 

(Pipo) := A{po/8, zi) + E{po/8, zi) + F(po/8, zi) < e' . 
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By using (I3.20p and (13.250 with p = pq/8, we obtain inductively that 

^(efpo) = Aie^po/S, zi) + Ei9fpo/8, zi) + FiOlpo/S, zi) < e' 

(holding for A; = 1,2, ...). It then similarly follows from ()3.20p and ()3.25p 
that 

iVo) < -Me\~'po) + iVi(^i'~Vo)', (3.26) 



4 

where we have used the estimate 

G(^^Vo/8,zi) + GV2((?^-Vo/8,zi)<iV(||/||L6(Q(.o,Po/2))(^t^Vo)'. 

Now we use a standard iteration argument to obtain the decay rate of 
We iterate p.26p to obtain 

where we have used that 9i < 3/4. Since p £ (0,po/32) we can find k such 
that e'l^ <4:p< e^-^f^. Then 



where N = iV((9i, (/>(po), iVi, Po) and ao = ^"^g^^/'^^^ > 0- This yields ([3:21]) 

for the terms j4, i?, and F. The inequality for C(p, zi) follows from (j3.18p 
and the inequality for D(p,zi) follows by (j3.4p . □ 

3.3. Step 3. In the final step, we are going to use a bootstrap argument to 
successively improve the decay estimate ()3.24p . However, as we will show 
below, the bootstrap argument itself only gives the decay of E{p) no more 
than p^^^, for instance, one can get an estimate like 

\Vu\'^dz < Np^+i 

Q{zi,p) 

for any p sufficiently small. Unfortunately, this decay estimate is not enough 
for the Holder regularity of u since the spatial dimension is four (so that we 
need the decay exponent 4 + e according to Morrey's lemma). We shall use 
parabolic regularity to fill in this gap. 

First we prove Theorem ll.2[ We begin with the bootstrap argument. We 
will choose an increasing sequence of real numbers {ofclfcLi ^ (^0,5/3]. 

Under the condition (j3.24p . we claim that the following estimate hold 
uniformly for all p > sufficiently small and zi £ Q(zq, po/8) over the range 

A{p, zi) + E{p, zi) < Np'^^ C{p, zi) + D{p, zi) < Arp3a./2_ (3 27) 

We prove this via iteration. The A; = case for (I3.27P was proved in ()3.24p 
with a possibly different exponent ao- 
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Now suppose (13 .241) holds with the exponent a^. We first estimate A{p, zi) 
and E{p,zi). Let p = 6p where 9 = pf^, p = p^~^ and p G (0,1) to be 
determined. We use (j3.17p . (I3.27j) (for a^), and (I3.13j) to obtain 

Choose p = , then ()3.27p is proved for A{p)+E{p) with the exponent 

10 + ak 

of 

Qfc+i := min ^2p + 0^(1 - /i), ^0^(1 ~ p) - Sp, 4(1 - p) - 
12 

-ak G {ak,2). 



10 + Ok 

Then the estimate in (j3.27p (with d^+i) for C{p, zi) follows from (j3.18p . To 
prove the estimate in (I3.27P (with Qfc+i) for D{p, zi) we will use Lemma [3.41 
From (|3.3p and (j3.13p . we have 

D{-fp, z,) < N[j-''pl^''+^ + 7'/'^(p, zi) + j-^p^] . 

For any r small, we take the supremum on both sides with respect to p G 
(0, r) and get 

sup L»(7p,zi) < iV7-Vi°'=+i +iV7^/2 sup D{p, zi) + j'^r^. 

pe(0,r] pe(0,r] 

Now set 

5 

afc+i = min{-,dfc+i}. (3.28) 

By using a well-known iteration argument, similar to ()3.26p . we obtain the 
estimate in (I3.27|) (with a^+i) for D{p). And the estimates in p.27p (with 
Ofc+i) for A{p) + E{p) and C{p) are still true. Then we have shown how to 
build the increasing sequence of {ofc} for which (I3.27[) holds. 
Moreover, 

10 , , 10 , 

2 - Ofc+i = 2 - ttk) < — — — 2 - ak). 

10 + afc 10 + ao 

5 

Thus, we can find a m that am = - according to (I3.28[) because otherwise 
afc = cifc — )• 2 as A: — )• oo. 

We have got the following estimates via the bootstrap argument: 



sup 



\u{x,t)\'^ dx < Np'^^^, (3.29) 



ti-p2<t<ti JB{xi,p) 

[ \uf + \p- [pUJ^^^ dz < Np^+I. (3.30) 
jQizup) 

Now we rewrite (11. ip (in the weak sense) into 

dtUi - Aui = -dj{uiUj) - dip + fi. (3.31) 
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Finally, we use the parabolic regularity theory to improve the decay esti- 
mate of mean oscillations of u and then complete the proof. Due to ()3.29p 
and (|3.30p , there exists pi E (p/2,p) such that 

'B(xi,pi) Jti-pl Js(xi,pi) 

(3.32) 

Let V be the unique weak solution to the heat equation 

dtv - Av = in Q{zi,pi) 

with the boundary condition Vi = Ui on dQ{zi, pi). It follows from the 
standard estimates for heat equation, Holder's inequality, and (|3.32p that 

sup |Vv| 

Q{zi,Pi/2) 

ftl [■ [■ 



<Np:^^'^ I I \v\dxdt + Npl^ I \v{x,ti - p{)\dx 

Jti^p'j Js{xupi) Jb{xupi) 



<Np-^+li. (3.33) 
Denote w = u — v. Then w satisfies the linear parabolic equation 

dtWi - Awi = -dj{uiUj) - diip- [p]xi,p) + fi 

with zero boundary condition. By the classical Lp estimate for parabolic 
equations, we have 

+ ^^/5i||/IIl3/2(Q{^i,Pi)), 

which together with (j3.30p and the condition / G L6,ioc yields 

f \Vwf/^dz<Np'^^/^. (3.34) 

Jq{zi,pi) 

Since \Vu\ < |Vw| + \S/v\, we combine p.33p and (I3.34D to obtain, for any 
r G (0,p/4), that 



/ \Vuf/^dz<Np''/' + 

jQ{zi,r) 



'Q{zi,r) 

Upon taking r = p^^/^'^/4 (with p small), we deduce 



[ \Vu\^/^dz < Nrf^, (3.35) 

jQ{zi,r) 



where 



n 132 ^ 3 
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Since u e ^3/2 loc ^ weak solution to ()3.3ip . it then follows from Lemma 
13.21 (j3.35p . and (j3.30p with r in place of p that 

jQ{zi,r) 

< iVr3/2 [ |Vn|3/2 + (|n|2)3/2 + |p _ [p],^^,|3/2 + r^/^\ff/^) dz 

JQ{zi,r) 

By Campanato's characterization of Holder continuous functions (see, for 
instance, \1Q\ Lemma 4.3]), u is Holder continuous in a neighborhood of zq. 
This completes the proof of Theorem 11.21 

Theorem 11.11 then follows from Theorem 11.21 by applying Proposition 13.61 
Finally, Theorem 11.31 is deduced from Theorem 11.11 by using standard argu- 
ment in the geometric measure theory, which is explained, for example, in 
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